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ABSTRACT 

Debris discs are known to exist around many planet-host stars, but no debris dust 
has been found so far in systems with transiting planets. Using publicly available cat- 
alogues, we searched for infrared excesses in such systems. In the recently published 
Wide-Field Infrared Survey Explorer ( WISE) catalogue, we found 52 stars with tran- 
siting planets. Two systems with one transiting "hot Jupiter" each, TrES-2 and XO-5, 
exhibit small excesses both at 12 /im and 22 /^m at a > Str level. Provided that one or 
both of these detections are real, the frequency of warm excesses in systems with tran- 
siting planets of 2-4% is comparable to that around solar-type stars probed at similar 
wavelengths with Spitzer^s MIPS and IRS instruments. Modelling suggests that the 
observed excesses would stem from dust rings with radii of several AU. The inferred 
amount of dust is close to the maximum expected theoretically from a coUisional cas- 
cade in asteroid belt analogues. If confirmed, the presence of debris discs in systems 
with transiting planets may put important constraints onto formation and migration 
scenarios of hot Jupiters. 

Key words: planetary systems - planets and satellites: detection - circumstellar 
matter ~ stars: individual: XO-5, HAT-P-5, TrES-2 parent star, CoRoT-8 



1 INTRODUCTION 

Many debris discs have been found in systems with known 



radial velocity (RV) planets (e.s. Beichman et al.l 2005a 



Moro-Martm et al 



20091 : iKospal et al" 



l2007l: iTrilling et al.. 2008; .Brvden et al.1 
20091 ) ■ and a few systems with debris 



discs and dir ectly ima ged plane tary candidates are known 
l|Kalas et al.i r2008: Mar ois et ai] |2008: LaEiango ct al. 20l3)- 
However, debris dust has not been found yet in systems with 
planets detected by transits. 

In this paper, we search for debris dust in systems 
with transiting planets, using publicly available catalogues 
of transiting planets and several infrared (IR) surveys. The 
motivation is obvious. A successful search would extend the 
list of known "full" planetary systems that harbour both 
planets and asteroid or Kuiper belt analogs. Furthermore, it 
is the transit technique that allows determination of many 
planetary parameters, such as masses, radii and densities, 
and can provide insights into properties of planetary atmo- 
spheres and interiors. Finally, transiting planets are on the 
average even closer to their parent stars than those discov- 
ered by the RV method, which might be related to somewhat 
different formation circumstances. Therefore, systems with 
transiting planets are of special interest. Detection of plan- 
etesimal belts, which are leftovers of planet formation, could 
help constraining various formation and evolution scenarios 



of those planets. And conversely, precise knowledge of plan- 
etary parameters could put constraints on the properties of 
the planetesimal belts. For instance, accurate masses and or- 
bits of planets would result in tighter constraints on dynam- 
ical stability zones and thus location of planetesimal belts. 



2 SEARCH FOR DUST 

A list of 93 currently known systems with tr ansiting 
plane ts was taken from exoplanets.org (Wright et al.l 
l201ll fl. This Ust was com pared with target lists of sev- 

Jl98i), 



ISO 



eral IR missions: IRAS (iNeugeb auer et al.l . _ _ . . , _ _ 
(Kc ssler et"aLlll996l), S pitzer l| Werner et al.ll2004l), AKAR I 
(.Murakami et al.l |2007| ). and WISE (* Wright ct al.' '2010l), 



which we accessed through IRSA, the NASA/IPAC Infrared 
Science Archive at |http : //irsa. ipac . caltech. edu| 

2.1 IRAS, ISO, Spitzer, and AKARI 

Nearly all of the transit planet host stars are located at hun- 
dreds of parsecs from the Sun and are thus faint. Accord- 
ingly, we had not expected to find them in older, and shal- 
lower, IRAS and ISO catalogues and indeed, have not found 
any. For example, of 93 systems with transiting planets listed 
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in [exoplanets . org , only five are within 50 pc. These are 
GJ 436, GJ 1214, HAT-P-11, HD 189733, and HD 209458. 
None of them appears in IRAS, ISO, and WISE catalogues. 
Three of them, GJ 436, HD 189733, and HD 209458, have 
been probed by Spitzer/MlPS at 24 a nd 70 ^m, yielding no 
excess detection (|Brvden et al.l |2009| ) . We found an entry 
for the latter star in the AKARI catalogue, reporting a de- 
tection at 9 ^m, which is consistent with the photospheric 
level. Note that HD 209458b was the first exoplanet found t o 
transit the disc of its parent star jCharbonneau et al.ll2000l ). 

We have also identified two more distant transit planet 
host stars that were observed by Spitzer/MIPS: HD 80606 
and HD 149 026. No excess at 2 4 and 70 /im was found for 
HD 189733 (|Brvden et al.ll2009t ). For HD 80606, the result 
is am biguous due to pointing problems (jCaroenter et al.l 
I2OO8I ). AKARIYias observed one more transit host star, too: 
HD 149026. It has been detected at 9 pm, showing no excess. 

2.2 WISE 

The search in the WISE Preliminary Source Catalog was 
more successful. This is perhaps not a surprise, given the 
broad sky coverage (57%) of the catalogue, an excellent sen- 
sitivity of the instrument and thus a huge number of sources 
observed (257 million). The W/SiS catalogue provides mea- 
sured magnitudes in four bands Wi (i = 1, . . . , 4), which 
are centred at 3.4, 4.6, 12, and 22 iim. Of 93 systems with 
transiting planet candidates listed in [exoplanets ■ org] we 
found 53 with entries in the WISE catalogue. One source — 
CoRoT-14 — is irretrievably contaminated by ghost images 
in bands Wz and 14^4 and was excluded from further analysis. 

To select possible IR excess candidates amongst the re- 
maining 52 sources, we first converted the observed magni- 
tudes in the four bands into spectral fiux densities. Since no 
excesses are expected in bands W\ and W2, we made simple 
photospheric predictions for W3 and Wa from the W\ and 
W2 fluxes. At first, we roughly corrected the Wi and W2 
fiuxes for an expected average level of interstellar extinction. 
Considering that systems in our sample are typically at a few 
hundreds parsecs, we set Av to 0.5™, which translates to 
MWi) = 0.029™ and A{W2) = 0.012™ (|Rieke fc Lebofskvl 
1 19851 ). We then fltted the corrected W\ and W2 fluxes with 
a power law Fphot ~ F'^^^^\~^ , with -Fphot and h being the 
fitting parameters. Subtracting the expected photospheric 
fiux from the observed one, we derived the "excess fiux" 
F = Fobs — fphot in bands W3 and W4. The net uncertainty 
of a photometric point for a given star in the band Ws or 
W4, was computed as cr = ^o-phot + '^obs + '^Li- Here, aphot 
is the photospheric uncertainty, which we estimated from 
the combined uncertainties of the measurements in bands 
Wi and W2, given in the WISE catalogue. Next, f7obs is 
the measurement uncertainty in the bands of interest, W3 
and W4, also taken from the WISE catalogue. Finally, (Tcai 
is the absolute calibration uncertainty of the WISE instru- 
ment (2.4, 2.8, 4.5, and 5.7% for bands from Wi to Wa). The 
significance of an excess can now be defined as x = F/a. 

The distributions of x-values in the sample are shown 
in Fig. [T]for bands W3 (top) and W4 (bottom). The W3 
histogram appears close to gaussian, without any obvious 
outliers. However, the W4 histogram uncovers a bin con- 
taining four systems with x > 1.75(j, clearly separated from 
the gaussian bulk. These are XO-5, HAT-P-5, TrES-2, and 
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Figure 1. Histogram of x = F/c values at 12 fim (top) and 
22 fim (bottom) for systems with transiting planets. For compar- 
ison, overplotted are closest gaussian distributions, which have 
variances of 0.65 (top) and 0.28 (bottom). That these variances 
are smaller than unity may suggest that the procedure of calcu- 
lating cr, described in the text, is too cautious and overestimates 
the actual noise. The rightmost bin in the W4 histogram shown 
in red contains four excess candidates. We checked that slight 
modifications to the photospheric extrapolation procedure or ex- 
tinction correction (i.e. using Ay values between 0.0™ and 1.0™) 
slightly skew and shift the W3 histogram, but preserve the W4 
one, including the bin with the four outliers. 



CoRoT-8. These four excess candidates will be checked in 
Sect.|3]more thoroughly, including an in-depth photospheric 
analysis and more accurate uncertainty estimates. 



3 ANALYSIS OF EXCESS CANDIDATES 

For the excess candidates, we have collected stellar data 
(Tab. [l| as well as optical and near-IR ph otometry. In 



the visual, we used the USNO-Bl.O Catalog ( Monet et al 
200a) . the Guide Star Catalog, Vers. 2.3.2 ( Lasker et al 
120081 ) ■ and the Afl-Sky Com piled Catalogue of 2.5 Mil 



lion Stars (jKharchenko fc Ro cscr 2009). The near-IR data 
stem from the 2MAS S AU-Sky Catalog of Point Sources 
(jSkrutskie et all I2OO6I ). For transforming the B, V, R, I 
magnitudes into units of fiux density [Jy] , we used the John- 
son calibratia^ system and for the 2MASS J, H, Ks bands 
the lCohen et al. (2003) calibrations. 

At the first step, this photometry was corrected for 
interstellar extinction. Since for distances considered here 
the latter is known to correla te with distance only weakly , 
we used colour indices from iKenvon fc HartmannI (|l995[ ) 
and the spectral type as given in Tab. [l] to derive the 
be st-flt Ay from mu l tiple colours and then Ax/Av ratios 
of iRieke fc Lebofskvl (| 19851 ) to compute extinction for the 
wavelengths of all photometry points. The derived Ay val- 
ues are within 0.1™ for non-CoRoT stars, but as large as 
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Table 1. Stellar parameters. 



Star 




d [pc] 


SpT 


[A/q] 


Teff[K] 


R* [Rq] 


Ref. 


d [pc] 




XO-5 


12.1 


260 ± 12 


late G 


0.88 ± 0.03 


5370 ± 70 


1.08 ±0.04 


[1] 


286 (+10%) 


5800 (+8%) 


HAT-P-5 


12.0 


340 ± 30 


early G 


1.16 ±0.06 


5960 ± 100 


1.17 ±0.05 


[2] 


340 (±0%) 


6400 (+7%) 


TrES-2 


11.4 


230 


GO V 


1.08 ±0.11 


5960 ± 100 




[3] 


219 (-5%) 


6400(+7%) 


CoRoT-8 


14.8 


380 ± 30 


Kl V 


0.88 ± 0.04 


5080 ± 80 


0.77 ± 0.02 


[4] 


304 (-20%) 


5400 (+6%) 



References: [1] IPal et al.1 l|2009f ). [2] iBakos et all l|2007l '). [3] lO'Donovan et al.1 l|2006l ). [4] iBorde et all l|2010l ). 
Note: The last two columns list best-fit values and their deviation from the starting values. 




Figure 2. Spectral energy distributions (SEDs) of four selected 
stars. Grey solid line: predicted extinction-corrected photosphere. 
Diamonds and open circles: visual and near-IR photometry data 
before and after correction for interstellar extinction, respectively. 
Filled circles: extinction-corrected WISE data (errors bars are 
""obs)- Black dashed line: modified BB model. 



1.3™ ± 0.2™ for CoRoT-8. The extinction in W3 and W4 
bands does not exceed 0.04™ (CoRoT-8). 

At the second step, we performed a minimum fitting 
of the extinction-c orrected stellar photos pheric fluxes by 
NextGen models (jHauschildt et all 1 1993). only to data 
points between 1 /jm and 5 fim. This is because those 
wavelengths are short enough not to expect any excess 
emission, but are long enough for interstellar extinction to 
be small. The interval from 1 fj,m to 5 fim includes three 
2MASS points (J, H, Ks) and two WISE points (3.4 
and 4.6 /im), all of which were given equal weights. Since 
the surface gravity of our stars (logg between 4.33 and 
4.61) and their metallicity ([Fe/H] from -0.2 to +0.31) 
deviate from the solar values only slightly (see Tab. [1] 
for references), we assumed a log g of 4.5 and the solar 
metallicity. Using Toff and d listed in Tab. [T] as starting 
values, we varied the temperature by ±400 K in 200 K steps 
and the distance to the star within ±30% in 5% steps to 
derive best-fit values of these two parameters. This yielded 
deviations from the starting values of up to 8% in Tctt and 
up to 20% in d (Tab. [T]). The results with the photometric 
points overplotted are shown in Fig. [2] Importantly, TrES-2 
and XO-5 and possibly also HAT-P-5 reveal small excesses 
in band W3 as well, which were not seen in Fig. [T] (top). 

We now come to a detailed analysis of the fluxes and 
their uncertainties. Denote the observed flux by J'obsi the 
extinction-corrected one by Fobs, the predicted extinction- 
corrected photospheric flux by Fphot, and the excess flux 
by F = Fobs — Fphot. As in Sect. [2j the net uncertainty 
of F for a given star in band W3 or W4 is computed as 
a — y^o^hff^^^"^rf)^"^^"^Ml• "^^^ measurement uncertainty 
(Jobs and the calibration uncertainty tJcai are included as de- 
scribed before. However, the photospheric uncertainty cTphot 
is now a by-product of the fltting procedure. It is dominated 
by a scatter in J, H, Ks, Wi, and W2 points (the error bars 
of the points themselves as well as the uncertainty of the 
extinction correction are much smaller). All the quantities 
above, and the resulting excess signiflcance x ~ F/a, are 
listed in Tab. [5] Nearly all excesses are at ~ 2a level, whereas 
usually a 3a excess is treated as a significant detection. 
However, in the cases of TrES-2, XO-5, and HAT-P-5, the 
excess is detected in two bands. The combined multi-band 
{W3 and W4) gaussian statistics suggests the significance 
level for these sources of 3.28, 3.23, and 2.82, respectively. 
This finally selects two systems, TrES-2 and XO-5, as 
> 3(j-significant and thus the best excess candidates. The 
binomial probability that one of these two detections is false 
is only 6.4%, and the probability that both are false is as low 
as 0.2%. The expected number of false detections at > 3a 
level is just 0.14; we detected two excesses at that level. 
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Table 2. Fluxes [mjy], uncertainties [mjy], and significance of excesses. 



System Band F*^^ Fobs -F'phot F o-phot c^ohs cTcai cr X Xjoint 



TrES-2 


W3 


3.71 


3.72 


3.10 


0.62 


0.20 


0.11 


0.17 


0.29 


2.17 


3.28 




W4 


2.22 


2.22 


0.93 


1.29 


0.12 


0.58 


0.13 


0.61 


2.12 


XO-5 


W3 


2.38 


2.38 


1.92 


0.47 


0.09 


0.14 


0.11 


0.20 


2.35 


3.23 




W4 


2.16 


2.16 


0.58 


1.58 


0.08 


0.85 


0.12 


0.86 


1.84 


HAT-P-5 


W3 


1.99 


2.00 


1.74 


0.26 


0.09 


0.10 


0.09 


0.16 


1.58 


2.82 




W4 


1.93 


1.94 


0.53 


1.41 


0.09 


0.68 


0.11 


0.70 


2.03 


CoRoT-8 


W4 


2.17 


2.24 


0.24 


2.00 


0.09 


0.96 


0.13 


0.97 


2.05 





Columns: Observed flux i^*]-,s' observed flux after correction for extinction -Fobsi expected photospheric flux i^photi excess flux _F; 
uncertainty of photospheric flux tTphot i observation uncertainty (Tobs i absolute calibration uncertainty Ccal i net uncertainty of the excess 
flux o", excess significance level in a single band x; joint (Wis and Wa) significance level Xjoint- 



Table 3. Dust parameters inferred from the observed excesses and parameters of transiting planets. 



System 


Blackbody 






Modified blackbody 






Planet 










■^blow 


so 






Md 


fd 






Mpi 




[K] 


[AU] 






[K] 


[AU] 






[AU] 




[Mj,p] 


TrES-2 


218 


1.7 


0.4 


2.1 


155 


5.8 


5 X IQ-^ 


3 X 10-"' 


0.037 


(fixed) 


1.28 


XO-5 


181 


2.2 


0.4 


2.0 


133 


8.0 


1 X 10-4 


6 X 10-4 


0.051 


0.049 


1.06 



Note: Planetary parameters are taken from the papers listed in Tab. [T] 



4 PRESUMED DUST BELTS 

In what follows, we estimate the parameters of dust that 
would produce the excesses in the best candidate systems, 
TrES-2 and XO-5, provided these are real. Since the ex- 
cesses are of low significance and the data are limited to 
two photometry points, a detailed SED modelling based on 
various assumptions about the size distribution and com- 
position of dust is not warranted. Instead, we used a pure 
blackbody (BB) and a modified BB emission model. In the 
latter case, we assumed a single grain size so and the opac- 
ity index of —2 beyond A = 27rso- The effective grain size so 
was chosen in the following way. Assuming a power-law size 
distribution with the index g = 3.5 and the lower cutoff ra- 
dius Smin of twjce the radiation pressure blowout limit Sbio w 
(see. e.g.. iKrivov et al.ll2006l : iThebault fc Augereaull2007^ ■ 
we have equated the emission of a disc of grains having such 
a size distribution and the emission of a disc composed of 
equal-sized grains of radius sq: 




abs 

(A,s) B^{\,Td{rd,8)) s^-" ds 



= Q!A>s{>^,so) B„{\,Td{rd,so)) si I s '' ds, (1) 

where rd is the distance from the star, A is the wavelength 
where excess emission is observed, is the Planck intensity, 
(3abs(A, s) is the grain absorption efficiency, and Td is the 
grain temperature. In calculating Sbiow, we assumed the unit 
radiation pressure efficiency and the bulk density of 3 g cm""^ 
and took the stellar parameters from Tab. [T] Equation ^ 
was solved for so. 

We then sought pure BB and modified BB curves that 
reproduce F(12jLtm) and F{22fj,m). This has yielded esti- 
mates of the temperature, location, mass, and the fractional 
luminosity of the emitting dust (Tab. [3]). When deriving the 
dust mass, we converted the mass of grains with s = sq 



into the mass of grains with s < 1 mm, assuming a power- 
law size distribution with a slope of 3.5. Both stars, TrES-2 
and XO-5, appear to have rings with radii of 6-8 AU and 
fractional luminosities in the range (3-6) X 10"*. We stress 
that all these inferred values are quite uncertain, because 
they rest on scarce photometric data and their derivation 
involves a number of simplifying assumptions and poorly 
known parameters. 



5 CONCLUSIONS AND DISCUSSION 

We have found that two out of 52 systems with transiting 
planets observed by WISE reveal warm two-band (12 and 
22 /xm) IR excesses at > 3a level. The probability that both 
excesses are real is 94% and that one of them is real is 99.8%. 

Provided that one of the two systems, or both, do pos- 
sess a warm disc, this would imply the excess incidence rate 
of 2-4%. For comparison, the Spiteer/MIPS detection rate 
of 24 nm excesses aro und old (^4 Gyr) fie ld stars was found 
to be 1/69 (1 ± 3%) (|Brvden et al.ll2006l ). Another sample 
of solar-type stars probed by MIPS at 24 ^ m resulted in 
fa 4% detection rate , avera ged over all ages (|Trilling et al.l 
aOOi ). iLawler et all (|2009l ) analyzed Spitzer /IRS observa- 
tions of nearby solar-type stars and found excess around 
12% of them in the long- wavelength IRS band (30-34 /im), 
but only 1% of the stars have detectable excess in the short 
wavelength band (8.5-12 fim). Thus the frequency of warm 
excesses around solar-type stars with transiting planets seen 
in the WISE data may be comparable to that in unbiased 
samples of similar stars found with Spitzer. 

Each of the two systems discussed here hosts one known 
close-in planet and, if the excesses are real, an asteroid belt- 
size dust ring well outside the planetary orbit. In both cases, 
more planets could orbit both inside and outside the belts. 
Additional planets at < 10 AU could be revealed by in-depth 
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RV analyses, by trans its, or by transit time var iations of al- 
ready known planets (jMacieiewski et al] 2011al'l . The latter 
method was used for TrE S-2 fRatz et all 120091 1 and XO-5 
l|Macieiewski et al]|2011bl ). Non-detection is consistent with 
the presence of debris belts at several AU, which are incom- 
patible with plan ets in that region. In the case of TrES-2, 
iRatz et alj (|2009l ') noticed a second dip in the light curve, 
both in their own light curves and those published in the 
literature. This second dip has been observed several times 
and then disappeared. In addition to o ther possible reasons 
for this effect, discussed bv lRatz et al.l . it could be due to an 
occultation by material in the debris disc. Estimates show 
that a clump of dust produced in a recent collision of two 
~ lOOkm-sized planetesimals would bear enough cross sec- 
tion to account for such a dip before it is azimuthally spread 
into a ring in a few years, although the probability of wit- 
nessing such an event is low. 

Planets at largest orbital radii (> 10 AU) will be hard 
to find by the transit technique. Direct imaging and astrom- 
etry are not feasible either, since these systems are too old 
and too distant. It will also be difficult to search for possible 
Kuiper belt analogues on the periphery of the systems, be- 
cause they are too faint for far-IR facilities such as Herschel. 
However, future mid-IR instruments such as JWST/MIRI 
should have enough sensitivity to study warm dust in great 
detail, including dust grain spectroscopy. They may also take 
a closer look at further potential excess candidates such as 
HAT-P-5 and CoRoT-8 identified in this study. 

The origin and the production mechanisms of the 
presumed dust are unclear. We have computed the dust 
mass expected to be produced through a steady-state coUi- 
sional cascade in a belt of "aste roids" with modera te ec- 
centricities, using the model of iLohne et al.l (|2008l ) with 
a velocity-dependent critical fragmentation energy from 
IStewart fc Leinhardtl ^09). At ages of ~ 1 Gyr, the maxi- 
mum expected dust mass is ~ lO~*A/0 at rd — 10 AU and 
~ 1O~^M0 at rd = 6AU. Comparing with Tab. [31 we con- 
clude that the amount of dust in our systems is close to, or 
even somewhat greater than, the theoretical maximum al- 
lowed by a steady-state coUisional cascade. This means that 
we might have a similar difficulty that exists in explaining 
other systems with hot excesses that have been known be- 
fore, such as HD 69830 (|Beichman et al.ll2005bl l. Proposed 
scenarios for such systems include: supply of comets from 
an outer massive cometary reservoir, possibly following a re- 
cent dynamical instability such as the Late Heavy Bombard- 
ment; the inward-scattering and desintegration of a large 
object from such an outer reservoir; a recent major colli- 
sion between two large planetesimals fsee lPavne et al.ll2009l , 
and references therein). Finally, a possibility of a steady- 
state coUisional dust production can be resuscitated if one 
allows the asteroid s in the belt to have very eccentric orbits 
IWvatt et al.ll2010l ). Such a belt could result from shepherd- 
ing and scattering of an initial plan etesimal belt durin g the 
inward migration of "hot Jupiters" ijPavne et alj|2009l ). 
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